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ABSTRACT 


This  thesis  contains  a  handout  covering  the  Fundamental 
Physical  Laws  (Continuity,  Momentum,  and  Energy)  used  in 
Aeronautical  Engineering  which  are  transformed  from  the  con- 
trol mass  or  system  form  into  the  control  volume  form.   It  is 
intended  that  this  handout  serve  as  a  self-studying  guide  for 
students  in  the  core  of  the  Aeronautical  Engineering  Program 
at  the  Naval  Postgraduate  School  and  as  a  reference  during 
the  graduate  level  courses. 
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I.   INTRODUCTION 


A.   THE  PROBLEM  UNDER  CONSIDERATION 

In  the  Aeronautical  Engineering  Curricula  at  the  Naval 
Postgraduate  School  the  development  of  expressions  for  the 
Fundamental  Physical  Laws  (Continuity,  Momentum,  and  Energy) 
is  covered  at  different  depths  in  courses  such  as  Basic  Fluid 
Dynamics,  Basic  Thermodynamics  and  Gas  Dynamics  during  the 
core. 

Due  to  the  lack  of  time,  these  topics  are  usually  treated 
in  an  introductory  fashion  and  the  notes  taken  by  the  students 
are  necessarily  incomplete  with  some  repetition  between  the 
various  courses. 

Also,  in  some  cases  the  development  of  these  equations  is 
done  without  a  systematical  recognition  that  control -vol ume 
approach  is  used  in  opposition  to  the  system  or  fixed-mass 
approach  (because  of  the  inconvenience  of  the  latter  for  prob- 
lems related  to  fluids).   In  addition  to  the  problem  previously 
described,  it  happens  sometimes  that  the  student  is  not  aware 
that  the  equation  the  professor  is  using  is  a  "disguised"  form 
of  one  of  the  Fundamental  Physical  Laws  for  which  he  may  al- 
ready have  some  intuition. 

The  situation  may  be  depicted  in  the  following  chart: 


TOPIC: 

FUNDAMENTAL 
PHYSICAL  LAWS 
INTEGRAL  FORMS 

FUNDAMENTAL 
PHYSICAL  LAWS 
DIFFERENTIAL 
FORMS 

PARTICULAR  FOKMS 

MAVIER-STOKES, 

ETC. 

DEVELOPED 
AND 
USED  IN 
COURSES: 

2041,  2042 

3043 

3043,  4431  ,  4432 

4501 ,  4502,  4511 

4512 

STAGE  IN 

THE 
PROGRAM: 

CORE 

GRADUATE 

OBJECTIVES: 

CREDIBILITY 

^ AND  ^  ^CONTINUITY  OF  SUBJECT  MATTER-^ 

RIGOR 

HflMnTMr'''"                   — •) 

-< '    n  M  IN  U  i  1  >!  L  o  o 

B.   TYPE  OF  STUDENTS  ATTENDING  THE  CORE 

A  considerable  number  of  students  taking  the  sequence  of 
core  courses  are  coming  from  the  fleet  after  a  number  of  years 
(between  four  and  eight),  since  the  completion  of  under- 
graduate education  in  Science  or  Engineering. 

It  can  be  said  that  sufficient  background  exists  in  most 
of  the  cases  but  it  is  buried  under  several  years  of  perfor- 
mance of  nonacademic  duties. 

The  most  noticeable  effect  of  these  inactive  years  is  the 
difficulty  in  dealing  with  abstract  concepts  such  as  mathe- 
matical symbols  and  models  intended  to  depict  a  situation  of 
engineering  interest  (the  control  volume,  for  example). 


C.   A  PROPOSED  SOLUTION  FROM  A  STUDENT  STANDPOINT 

The  author  believes  that  a  handout  received  by  the  students 
before  or  during  the  core  will  provide  them  with  a  set  of  notes 


for  self-study  during  the  core  and  for  future  r-3-^^rence  in 
the  graduate  sequence  of  courses. 

The  proposed  handout  on  the  Fundamental  Physical  Laws  is 
a  version  of  material  covered  in  many  books,  some  of  them 
being  more  rigorous  and  elegant  but  often  not  appealing  to 
the  students. 

What  is  it  that  bothers  a  student  in  this  situation  when 
the  topic  is  treated  in  class  or  in  his  textbooks? 

First  of  all,  the  terminology.   The  reader  must  be  warned 
that  the  author  is  perhaps  strongly  biased  in  this  respect 
because  English  is  not  his  native  language.   But  even  students 
without  this  problem  sometimes  are  confused  by  words  or  do  not 
know  the  exact  meaning  of  some  terms  in  Thermodynamics  or 
Fluid  Dynamics.   For  example,  look  up  the  word  "steady"  in  a 
dictionary.   Reference  1  gives  the  following  meanings  for  it: 
firm,  stable,  not  shaky,  regular,  uniform,  continuous,  constant 
in  behavior,  etc.   Assuming  that  the  student  knows  all  these 
meanings,  if  the  word  "steady"  is  used  without  rigorously 
defining  it  in  a  scientific  context  (i.e.,  partial  derivatives 
with  respect  to  time  are  zero),  a  sensation  of  ambiguity  will 
remain  and  no  convenient  abstraction  can  be  expected  to  last 
in  the  student's  mind. 

Important  terms  must  be  carefully  defined  from  the  begin- 
ning and  the  other  terms  introduced  as  needed  during  the 
development.   The  use  of  these  terms  must  be  shown  by  means 
of  examples  and  the  degree  of  retention  by  the  student  must 
be  frequently  tested.   (Notice  how  the  term  system  is  treated 
in  the  text . ) 
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A  second  source  of  trouble  is  the  prolific  contribution 
to  Science  by  men  such  as  Newton,  Euler,  Bernoulli  (three  of 
them!   Jacob  (1664-1705);  Johan  (1667-1748);  and  Daniel 
(1700-1782));  Laplace  and  many  others.   When  someone  quotes 
"Newton's  Law,"  to  which  one  is  he  referring  —  the  First, 
Second,  or  Third?   Maybe  it  is  Newton's  Law  of  Friction. 
Sometimes  the  student  is  taking,  simultaneously,  a  course  in 
Mathematics  where  the  same  names  are  also  quoted. 

The  fundamental  points  for  a  student  in  the  core  in  this 
respect  are: 

1.  Not  to  get  confused  by  the  terminology  and  partic- 
ul ar  names . 

2.  To  recognize  where  an  equation  originates;  what  is 
the  significance  of  the  various  terms;  and  what  are  the  con- 
straints and  assumptions  for  its  derivation. 

3.  To  make  actual  use  of  the  equations. 

Another  problem  arises  from  Algebra.   It  is  not  that  this 
branch  of  Mathematics  is  obscure  or  difficult  for  a  student  at 
this  stage  of  his  education,  but  that  it  is  simply  uninteresting 

Oftentimes  during  a  long  and  tedious  derivation  of  an 
equation  the  student's  mind  wanders  and  the  sensation  of 
"magic"  results  when  a  final  and  neat  expression  appears.   The 
following  corrective  actions  may  be  proposed: 

1.  State  from  the  beginning  the  starting  and  ending 
points  of  the  derivation. 

2.  Underline  the  milestones  in  the  derivation. 

3.  Ask  questions  of  the  student  in  the  middle  of  the 
derivation  to  keep  him  interested  or  simply  awake. 


4.  Skip  on  purpose  some  steps  but  ask  the  students 
to  go  through  them  for  completeness. 

5.  Use  the  most  compact  notation  available  but  only 
after  being  sure  that  the  notation  is  properly  introduced  and 
that  this  greater  degree  of  abstraction  is  consistent  with  the 
level  of  the  student. 

Perhaps  many  of  these  problems  are  avoided  by  using  a 
textbook  writing  technique  called  "Programmed  Instruction." 
The  opinion  of  the  author  is  that  writing  a  handout  fully 
using  this  technique  is  out  of  the  scope  of  a  Master's  Degree 
Thesis,  but  some  of  the  features  of  this  technique  can  be 
incorporated  and  the  possibilities  of  improvement  in  this 
respect  are  unlimited. 

D.   TyE  VARIABLE  OBJECTIVES  OF  THE  HANDOUT 

The  objective  of  the  proposed  handout  is  not  necessarily 
to  teach  a  technique  or  a  skill  (such  as  the  use  of  the  Laplace 
Transform  to  solve  certain  differential  equations).   The  objec- 
tives of  this  handout  are  variable  from  the  student  standpoint, 
accordingly  to  what  stage  of  his  education  he  is  at. 

At  the  beginning,  during  the  core,  the  main  text  of  the 
handout  is  intended  to  satisfy  the  objective  of  CREDIBILITY 
and  RIGOR.   The  student  should  be  familiar  (and  if  he  has  for- 
gotten, the  wery    first  part  of  the  handout  will  serve  as  a 
refresher)  with  the  system  or  fixed-mass  approach.   The  hand- 
out will  lead  him  neatly  and  rigorously  to  the  control  volume 
approach  through  a  transformation  or  extension  of  the  system 
approach  . 
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This  needs  to  be  done  only  once  and  it  can  be  done  using 
a  generalized  notation.   Also  it  is  believed  that  a  self- 
study  text  has  the  potential  of  removing  the  usual  class- 
room prob 1  ems  . 

In  the  next  steps  the  student  is  led  through  specializa- 
tions of  the  control  volume  formulation.   The  students  who 
are  in  the  core  are  not  expected  to  carry  out  these  further 
specializations,  but  he  is  expected  to  do  so  in  his  graduate 
courses.   Here,  the  objective  of  the  handout  is  to  provide 
CONTINUITY  OF  SUBJECT  MATTER  through  the  whole  sequence  of 
courses. 

The  same  objective  is  pursued  in  the  extensions  leading 
to  the  differential  form  of  the  Fundamental  Physical  Laws  and 
to  other  particular  forms.   He  is  not  expected  to  carry  out  -^ 
these  derivations,  but  he  is  expected  to  understand  where  the 
equations  come  from  and  (eventually)  the  use  of  them  in  the 
context  of  the  graduate  level  courses.   Here  the  "mystery"  of 
the  equation  written  down  without  proof  might  be  removed  and 
the  student  uneasiness  with  an  advanced  subject  alleviated. 

The  entire  handout  should  serve  properly  as  a  review  and 
reference  during  all  the  stages  of  his  career  at  the  Naval 
Postgraduate  School.   The  objective  in  this  respect  is  HANDI- 
NESS.   This  is  perhaps  the  most  ambitious  aim  of  this  work. 


E.   THE  LEARNING  PROCESS 

This  thesis  has  been  written  according  to  the  following 
basic  ideas  with  respect  to  the  learning  process: 

1.   The  assumption  that  the  student  is  motivated;  he  has 
the  desire  to  learn. 


2.  Mo  mathematical  background  is  assumed  beyond  undergraduate 
level . 

3.  Certain  basic  vocabulary  is  assumed  but  the  important 

concepts  are  defined  or  explained  carefully  and  their 
understanding  emphasized  and  tested  in  several  parts  of 
the  text. 

4.  There  would  be  systematic  reading  of  the  context  during 
the  self-study  use  of  these  notes.   The  student  will  not 
proceed  to  another  section  without  a  fair  understanding  of 
the  previous.   (Approximate  relative  times  for  reading 

are  stated . ) 

Along  the  text,  several  quizzes,  programmed  questions  and 
exercises  are  inserted  to  help  the  student  recognize  his  own 
degree  of  retention  of  what  he  was  supposed  to  learn. 

When  a  nev^  quantity  is  introduced  and  assigned  a  symbol, 
the  corresponding  units  are  shown  in  the  English  system  only 
as  an  example  to  help  the  reader  recognize  and  associate  the 
nature  of  the  quantity. 

Also,  to  some  restricted  extent,  some  provisions  are  made 
in  order  to  allow  the  student  to  extend  his  recently  obtained 
knowledge  to  new  situations. 

F.   SUMMARY  OF  THE  PRESENTATION  OF  THE  SUBJECT  IN  THE  HANDOUT 

The  general  idea  followed  during  the  development  of  the 
thesis  was  to  treat  the  topic  at  the  highest  degree  of 
generality  and  compactness  consistent  with  the  level  of  the 
student  and  then  to  particularize  to  arrive  at  useful  forms 
of  the  equations. 
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It  can  be  mentioned  in  this  respect  that  the  reduction 
of  the  expressions  of  the  three  Fundamental  Physical  Lav/s  to 
one  using  a  generalized  notation,  was  based  on  Ref.  2  after 
consulting  several  undergraduate  level  books  and  notes 
[Refs.  3  through  11]. 

The  concepts  involved  in  the  control  volume  approach  were 
applied  to  this  generalized  expression  of  the  Fundamental  Phys- 
ical Laws  for  a  system.   Later  a  general  expression  of  the 
Fundamental  Physical  Laws  for  a  control  volume  was  developed. 
Then  this  expression  was  specialized  to  the  particular  laws  in 
Integral  form  and  later  in  differential  form. 

The  advantage  of  this  method  is  that  the  mathematical 
transformation  needs  only  to  be  done  once,  enabling  the  student 
to  concentrate  on  the  actual  mathematics  without  having  to 
identify  the  "proof"  with  a  single  Fundamental  Physical  Law. 

The  steps  in  parti cul arizati on  of  the  equations  are  con- 
tained in  the  supplements  where  the  Navier-Stokes ,  Bernoulli 
and  Euler  equations  are  developed  having  in  mind  that  yery 
often  a  fair  knowledge  of  these  equations  is  required  in 
courses  such  as  Boundary  Layer  Theory,  Convective  Heat  and 
Mass  Transfer,  Magnetohydrodynami cs  ,  Advanced  Gas  Dynamics, 
Turbomachi nery  and  Aerodynamics  of  Wings  and  Bodies. 

G.   THE  SCOPE  OF  THE  APPROACH  CHOSEN 

This  thesis  is  a  student's  point  of  view  on  the  subject, 
taking  into  consideration  the  difficulties  found  during  his 
education . 
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This  type  of  approach  may  be  helpful  to  fellow  students 
in  the  same  v/ay  as  it  has  been  for  the  author. 

Under  the  title  of  this  thesis,  a  broader  range  of  sub- 
jects can  be  covered.   However,  some  amount  of  limitation  was 
found  to  be  necessary  in  order  to  fulfill  most  of  the  purposes 
enunciated  earlier  in  this  Introduction. 

It  is  necessary  to  mention  that  the  treatment  of  the 
Equation  of  State  which  is  complementary  to  the  Fundamental 
Physical  Laws  has  been  omitted.   The  reason  for  this  is  that 
this  topic  is  covered  in  a  thesis  written  at  the  same  time 
by  another  student  [Ref.  12]. 

II.   THE  HANDOUT 

In  order  to  have  the  proposed  handout  in  a  style  and 
format  more  useful  for  its  purpose,  it  is  inserted  in  the 
Appendix. 

Ill .   A  SURVEY  OF  THE  USE  OF  THE  HANDOUT  BY  STUDENTS 

An  earlier  version  of  the  handout  was  distributed  to  test 
its  acceptance  by  the  students. 

This  was  done  in  Course  AE  2041,  Basic  Thermodynamics, 
where  Ref.  6  was  used  as  a  textbook.   Upon  the  completion  of 
the  Course  the  instructor  made  a  survey  of  the  opinion  of  the 
students  about  the  handout.   The  following  are  some  results 
of  such  a  pol  1  : 
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95  %    read  the  handout  in  its  entirety. 

88  %    had  sufficient  background  to  understand  the 
material . 

95  %    felt  that  the  material  was  of  sufficient 
difficulty  to  merit  a  handout. 

95  %    felt  that  the  material  was  presented  in  a 
logi  cal  order. 

When  asked  whether  they  would  prefer  their  textbook,  the 
handout,  neither,  or  both,  the  response  was: 
40  %    present  text 
35  %    handout 
15  %    both 
10  %    neither 

One  important  conclusion  based  on  the  survey  is  that  at 
least  the  handout  received  was  a  desired  instructional  aid. 

With  the  improvements  made  since  then  and  other  improve- 
ments that  can  be  made  after  testing  it  again,  the  handout  may 
be  transformed  in  a  very  useful  instructional  aid  to  be  re- 
ceived by  the  students  at  the  beginning  of  their  career  at 
the  school . 
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I V .   RECOMMENDATIONS  FOR  FUTURE  WORK 

The  handout  may  be  improved  mainly  in  the  supplements. 

In  Supplement  C,  the  differential  form  of  the  Energy  Equa- 
tion is  only  stated,  suggesting  that  its  development  is  simi- 
lar to  the  development  of  the  Navier-Stokes  Equations  which 
in  turn  are  not  derived  in  detail  for  the  sake  of  conciseness. 

Perhaps  a  few  more  steps  may  be  added  in  the  development 
of  the  Navier-Stokes  Equations.   The  complete  derivation  of 
the  Energy  Equation  in  differential  form  is  needed  to  satisfy 
the  requirement  of  completeness  for  reference  use. 

There  are  also  opportunities  for  improvement  in  the  main 
body  in  the  style  of  the  presentation.   After  being  used  by 
the  students,  the  programmed  characteristics  may  be  increased 
to  overcome  difficult  steps. 
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I.   PRELIMINARIES 

A.   A  SUGGESTION  ON  HOW  TO  USE  THIS  HANDOUT 

This  handout  is  supposed  to  promote  a  better  understanding 
of  the  Fundamental  Physical  Laws  (Continuity,  Momentum,  and  En- 
ergy) and  is  to  be  read  out  of  the  classroom  (i.e.,  self-study). 

To  obtain  better  results,  (usually)  at  the  beginning  of 
each  chapter  or  section  there  is  a  note  suggesting  how  far 
you,  the  student,  should  go  on  reading  continuously.   We  do 
not  forecast  how  long  it  will  take  to  read  and  understand  each 
unit  because  of  individual  differences,  but  our  guess  is  that 
the  first  two  units  will  take  you  about  one  hour.   The  esti- 
mates on  the  following  sections  are  indicated  relative  to 
the  first. 

If  you  do  not  have  the  time  now  to  read  continuously  as 
suggested,  try  it  some  other  time,  because  partial  reading 
will  not  benefit  you  in  any  way. 

Also,  there  are  three  types  of  testing  material  inserted 
along  the  text  and  you  should  understand  their  purposes: 

1 .  Sel f-Che  ck  Qui  zzes 

The  purpose  of  these  quizzes  is  to  help  you  recognize 
your  own  degree  of  retention  before  proceeding  to  the  next 
page  of  the  material.   Be  sure  to  correct  the  wrong  answers, 
and  do  not  continue  unless  you  feel  you  understand  the  dif- 
ficulties you  had.   The  answer  is  usually  at  the  bottom  of 
the  page,  written  upside  down. 

2.  Programmed  Questions 

The  primary  purpose  of  these  questions  is  to  remind 
you  of  important  points  presented  in  previous  sections.  You 

must  fill  the  blank  space  left  on  purpose  and  the  answer  is 
usually  at  the  bottom  of  the  page,  written  upside  down. 
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3 .   Exercises 

These  dre    not  truly  testing  material.   ^heir  nurpose 
is  to  make  you  think  about  further  implications  of  the  subject 
in  an  open-ended  fashion.   No  anS'.yer  is  given  for  these. 

B.  FUNDAMENTAL  PHYSICAL  LAWS 

In  the  solution  of  problems  in  Aeronautical  Engineering, 
concepts  from  the  Physical  Sciences  called  "Fundamental  Phys- 
ical Laws,"  must  be  applied.   These  laws  relate  to  mass, 
momentum,  and  energy  and  are  commonly  referred  to  as  "conser- 
vation laws"  because  under  many  circumstances  they  are  truly 
conserved;  however,  the  word  conservation  in  its  common  usage 
takes  a  broader  meaning  which  may  be  stated  as  "accounted  for" 
or  "not  destroyed."   We  will  use  the  term  "Fundamental  Physical 
Laws  (FUPLAs)"  because  this  will  permit  us  to  include  the  second 
law  of  thermodynamics  relating  to  the  entropy  which  is  not  con- 
served in  a  natural  process.   In  addition  to  the  FUPLAs,  in- 
formation is  needed  about  the  fluid  itself  usually  contained  in 
the  Equation  of  State.   This  equation  describes  a  set  of  unique 
relationships  between  state  properties  and  may  be  found  in 
algebraic  form  (the  ideal  gas  law),  tabular  form  (the  steam 
tables),  or  graphical  form  (p-v-T  surfaces,  etc.). 

C.  BASIC  DEFINITIONS 

Read  continuously  from  here  to  the  end  of  Section  D. 

1 .   A  Word  about  the  Notation 

Through  all  this  work  a  ^    (caret)  will  be  used  over 
the  first  letter  of  a  word  which  has  a  specific  meaning  in 
this  text.   For  example,  whenever  the  word  cycle  appears  with 
a  ^    over  the  letter  "c"  (i.e.,  cycle),  means  a  precise  concept 
and  that  there  is  no  other  meaning  for  the  same  word. 

All  words  beginning  with  a  letter  accented  with  ^    are 
carefully  defined  in  the  following  list. 
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Fluid:   a  substance  that  deforms  continuously  under  the  action 
of  a  shear  stress. 

System:   a  collection  of  matter  of  fixed  identify  (Control 
Mass  )  . 

Control  Volume  (CV):   a  region  in  space  of  fixed  size  through 
which  fluid  flows.   It  may  be  stationary,  moving  uniformly 
(inertial),  or  accelerating  ( non- i nerti a  1 ) .   In  this  work  we 
will  consider  the  CV  always  fixed  with  respect  to  an  inertial 
coordinate  axis  system,   unless  otherwise  specified.   This 
volume  is  bounded  by  a  control  sur-^"ace. 

We  want  to  insist  on  the  idea  that,  although  the  con- 
trol volume  has  to  be  fixed  with  respect  to  an  inertial  sys- 
tem of  coordinates,  it  does  not  necessarily  need  to  be  fixed 
in  space.   The  following  illustration  shows  this.   The  vehicle 
may  be  moving  through  space  but  the  control  volume  is  fixed 
with  respect  to  the  vehicle. 


Figure  1.   Fixed  Control  Volume  with  respect  to  inertial 
frame  of  reference . 


Our  development  is  constrained  to  this  case  for  simplic- 
ity and  to  fit  into  the  scope  of  the  handout  like  this.   Impor- 
tant Aeronautical  Engineering  applications  require  non-inertial 
frames  of  reference,  such  as  in  turbomachi nery ,  for  example, 
[Ref,  13]. 
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Surroundings:   everything  outside  the  system 

Extensive  Property:   a  property  that  is  mass  dependent.   It' 
can  be  a  scalar,  as  energy  or  a  vector,  as  momentum 

Fundamental  Physical  Laws  (FUPLAs):   there  are  three  funda- 
mental physical  laws  which,  with  the  exception  of  rela- 
tivist ic  and  nuclear  phenomena,  apply  to  each  and  every  flow, 
independently  of  the  nature  of  the  fluid  under  consideration 

State:   the  condition  of  the  system  characterized  by  the 
values  of  its  properties 

Process:   the  path  of  the  succession  of  states  through  which 
the  system  passes 

Cycle:   a  system  that  undergoes  a  series  of  processes  and 
always  returns  to  its  initial  state  is  said  to  have  gone 
through  a  cycle 

Properties:   characteristics  of  a  system  that  define  its  state. 
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TABLE  I 


FUPLAs  (From  Observation) 

MATHEMATICAL  FORMULATION 

Conservat i  on  of  Mass 
Newton's  Second  Law  of  Motion 
First  Law  of  Thermodynamics 

Continuity  Equation 
Momentum  Equation 
Energy  Equation 

Note  that  Table  I  is  incomplete.   Other  FUPLAs  will  be  seen 
later. 

2.   A  Word  about  FUPLAs 


Perhaps  the  reader  has  become  uneasy  about  the  previous 
wording.   All  those  terms,  with  the  caret  (which  denotes  a 
specific  meaning)  and  the  solemn  title  "Fundamental  Physical 
Laws,"  may  lead  one  to  think  that  these  laws  have  something 
mysterious  or  difficult  to  understand. 

This  should  not  be  so. 

These  laws  were  deduced  from  observations  of  nature. 
Nature  behaves  this  way.   It  has  always  behaved  so  and  nobody 
has  reliably  observed  a  contradiction  to  this  behavior. 

This  is  the  only  "proof"  of  these  laws. 

Later  these  Fundamental  Physical  Laws  are  put  into  a 
mathematical  form  in  order  to  manipulate  and  work  with  these 
concepts  in  the  precise  and  synthetic  language  of  the  mathe- 
mati  cal  symbols . 

This  almost  trivial  example  shows  that  the  FUPLAs 
belong  to  everyday  experience:   If  you  bad  a  dozen  apples 
which  you  were  to  squeeze  to  make  apple  sauce  with,  you  would 
not  be  conserving  the  number  of  apples.   What  is  to  be 
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accounted  for  is  the  matter  (mass)  v;hen  the  seeds  and  other 
discards  are  taken  into  account. 

FUPLAs  are  the  generalization  of  many  reliable  observa- 
tions (i.e.,  matter  is  not  destroyed,  it  is  transformed  into 
other  forms  of  matter  or  into  energy). 
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SELF-CHECK  QUIZ  NO.  1 
Select  the  proper  answer 

1.  The  .Control  Volume  is: 

a.  a  collection  of  mass  of  fixed  identity 

b.  always  the  volume  occupied  by  the  system 

c.  a  region  in  space  of  fixed  size  through  which 
fluid  flows 

d.  always  stationary 

2.  The  Fundamental  Physical  Laws: 

a.  are  deduced  from  observation  of  nature 

b.  can  the  proven  using  highly  complicated  mathemat 
ical  procedures  and  as  a  result  of  this,  most  of 
the  textbooks  do  not  have  them 

c.  apply  to  certain  types  of  fluids,  flowing  under 
special  conditions 

d.  have  no  mathematical  expressions 

(e)2  pue  (D)L  :sj8msuv 


EXERCISES: 


What  is  conserved  when  an  object  of  mass  m  falls  from  a 
height  h  in  a  gravity  field? 

A  typical  reference  state  for  air  is  standard  conditions 
Do  you  know  what  this  means? 
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II  .   DEVELOPr-iENT  OF  EXPRESSIONS  FOR  THE  FUNDAMENTAL  PHYSICAL 

LAMS  FOR  A  SYSTEM 

Read  continuously  up  to  the  end  of  Section  E.  = 

A.   CONSERVATION  OF  MASS 

A  system  is  

Let  us  consider  the  following  system 


.—•  •  — •  •  —  •  ■—  •     •  •  • 

iv.;';.  GAS.  •••'-•ji 

' L  —  —  —  *—  —  Sly 

/r  — 

-■—■—■ 

System 


Figure  2.   Gas  in  a  closed  cylinder 

Experience  tells  us  that  matter  will  not  flow  in  or  out 
through  the  solid  boundaries  of  this  system  (or  any  other), 
as  time  elapses  or  as  the  identified  system  wanders  through 
space.   Also,  matter  will  not  be  created  spontaneously  within 
the  system.   From  this  last  statement  we  may  write 


m 


syst 


constant 


(1) 


^yst  -   ^syst^^  C2)- 

^:inu9pL  paxLj.  J.0  sseui  ^o  A:iL:^uenb  e   ijawsuy 
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If   p    =    ^   ,    then    m^^^^^    =      f^^^^odV 


(3) 


Differentiating    with    respect    to    time    in    Eq.    3    and    keeping 
in    mind    that  'the    left-hand    side    is    a    constant    (from    Eq.    1) 
we    get 


d_ 
dt  ^ 


syst 


pdV 


0 


(4) 


Note  that  in  the  above,  as  well  as  throughout  the  rest  of 
this  Section  II,  integrals  over  the  mass  (as  in  Eq.  2)  are 
being  used.   They  are  changed  into  volume  integrals  provided 
we  can  use  the  mass  density  p.   In  other  words,  for  the  system 
analysis,  integrals  over  a  fixed  quantity  of  mas  s  are  being 
dealt  with,  and  eventually  these  basic  laws  in  terms  of  control 
volume  will  be  expressed.   In  order  to  avoid  confusion  as  to 
which  one  the  volume  integral  belongs,  the  term  "syst"  will 
appear  by  the  integral  sign  when  a  fixed  mass  is  being  analyzed 
and  the  term  "CV"  when  a  fixed  or  control  volume  is  being 
analyzed . 

0 

B.   MOMENTUM  EQUATION 

Consider  a  differential  volume  element  of  our  system  which 
is  traveling  at  velocity  v.   This  volume  element  will  have  a 
mass  dm  and  agai n 


dm 
and  a  momentum 
d^ 
dP 


pdV 


vdm 


vpdV 


So  that  the  total  momentum  for  the  system  will  be 


syst 


d^ 


J 


J 


syst 


vpdV 


(5) 

(6) 
(7) 

(8) 


Differentiating  Eq.  8  with  respect  to  time 
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r 


r 


d_ 
dt 


> 


syst 


dP 


d_ 
dt 


J 


syst 


vpdV 


(9) 


The  Momentum  equation  is  the  expression  of  Newton's  Second 

Law  of  Motion  which  states: 

"The  time  rate  of  change  of  momentum  of  a  system 
is  proportional  to  the  net  external  force  acting  on 
system  and  takes  place  in  the  direction  of  the  net 
force .  " 

(Is  there  a  "proof"  for  this  law?) 


Therefore , 


r 


net 


_d 
dt^ 


dP 


(10) 


syst 

The  net  force  is  made  up  of  contributions  of  two  kinds  of  force, 
namely, 

(n) 


F   ..   =  >  F   T     +  >  F 

net    /       volume  /      surt 


ace 


The  above  distinction  regarding  the  net  force  is  needed 

because  of  the  different  nature  between,  say,  the  force  of 

gravity  and  a  pressure  force  (Supplement  A  expands  on  the  nature 
of  these  forces ) . 


Inserting  Eq.  9  into  Eq.  10: 


r 


net 


d_ 
dt 


J 


syst 


pv  dV 


If  there  is  no  net  external  force,  Eq.  8  becomes 


or 


d_ 

dt 


P 


r 


syst 


dP 


(12) 


(13) 


constant  (i.e.,  momentum  is  conserved) 

(14) 


The  student  should  eventually  learn  to  recognize  the  dif- 
ference between  the  collision  of  two  bill.iard  balls  and  of  two 
balls  of  putty  in  terms  of  the  momentum  equation. 
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Again,  all  that  has  been  considered  so  far,  v-n'th  respect 
to  the  Conservation  of  Momentum,  is  restricted  to  the  system 
approach.   At  this  point  the  above  will  not  be  'jery    useful 
if  the  system  is  difficult  to  identify. 

C.   ENERGY  EQUATION 

The  Energy  Equation  is  the  mathematical  expression  of  the 
First  Law  of  Thermodynamics,  which  can  be  stated  as  follows: 
"If  a  system  is  carried  through  a  cycle,  the  total  heat  added 
to  the  system  from  its  surroundings  is  proportional  to  the 
work  done  by  the  system  on  its  surroundings." 

Let  us  have  a  system,  that  is,  a  quantity  of  matter  of  fixed 
identity,  traveling  through  space  (changing  its  position  with 
respect  to  an  inertial  frame  of  reference  as  time  elapses). 

Let  the  path  of  the  system  be  in  such  a  way  that  it  will  re- 
turn  to  its  initial  state.   In  doing  so,  it  has  performed  a 
cycl e . 

During  the  cycle,  heat  may  be  added  to  the  system  by  its 
surroundings  and  in  turn,  work  may  be  done  by  the  system  on  the 
surroundings.   A  useful  cycle  is  depicted  in  Figure  3. 

Notice  that  the  previous  discussion  and  Figure  3  refers  to 

a  system  undergoing  a  cycle,  that  is,  

.   A  similar  discussion  may  be  made  for  a 


system  undergoing  a  single  process. 

Since  energy  is  related  to  heat  and  work,  for  a  process  — 
an  expression  for  the  First  Law  of  Thermodynamics  in  the 
rate  form   can  be  written: 


Rate  of  heat 
flow  into  the 
system 


rate  of  change  of 
energy  within  the 
system 


rate  of  work 
done  by  system 
to  the 
surroundings 

Here,  the  following  sign  convention  is  used:   Heat  added  to 

the  system  by  the  surroundings,  positive  (+).   Work  done  on 

the  surroundings  by  the  system,  positive  (+). 


(15) 


•8:tG:;s    [euL5LJ0    s:n.    o:^    Lua:^s/Cs    9L|:^ 
J.0    8:^e:^s   bi\%    suun^au    sApml^    ^^M^    sasssooud   j.o    saiuas    e 
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t> 
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out 


Figure    3.      Gas    cycle 


WORK 
W   (  +  ) 


Figure    4.      Sign    convention    for    the   First    Law 
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This  sign  convention  is  not  universal.   Some  authors  use 
other  conventions. 

Let  us  define  the  following  symbols: 

Q    rate  of  heat  flow  (e.g.,  in  ^^) 

^      sec 

IP   rate  of  work  done  (e.g.,  power  i.e.,  in  — ^—r-   or  wattsj 


sec 


Btu 


e   total  energy  per  unit  mass  (e.g.,  yf— ) 

Then,  one  may  write  the  energy  equation  for  a  system  from  Eq 
15:  r 


Q  -  IP  =  5z- 


dt_J  syst 


IP 


e  dm 


e  p  dV 


(16) 
(17) 


The  idea  that  energy  cannot  be  created  or  destroyed  within 
a  system,  but  that  it  can  be  exchanged  is  an  old  one  and  is 
easily  accepted  due  to  the  fact  that  it  has  never  been  ob- 
served to  be  violated  in  nature. 

The  problem  again  is  to  extend  the  Energy  Equation  to  cases 
where  it  is  difficult  to  follow  the  control  mass. 

D.   SUMMARY 

The  expressions  for  the  FUPLAs  for  a  system  is  summarized 
in  Tables  II  and  III. 

Table  II  will  be  extended  to  other  Laws  which  has  not  yet 
been  discussed. 

In  Table  III  there  are  some  symbols  not  defined  so  far. 
Refer  to  "Extension  of  the  Control  Volume  Approach  to  other 
Fundamental  Physical  Laws"  for  explanation  (Chapter  IV). 
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TABLE  II 


Fundamental 
Phys  i  cal 
Lav; 


Name  Usually  Given 
to  the  Equation 


Expression  for 
a  System 


Conservat i  on 
of  Mass 


Newton ' s 
Second  Law 
of  Motion 


Continuity  Equation 


d_ 
dt 


y 


^  pdV   =   0 
syst  ^ 


Momentum  Equation 

or 
Momentum  Theorem 


dt   syst 


vpdV 


net 


First  Law 
of 
Thermodynami  cs 


r 


Energy  Equation 


d_ 
dt 


y 


syst 


epdV 


Q  -  P 


TABLE  III 


Fundamental 

0 

Physi  cal 
Law 

Name  Usually  Gi  ven 
to  the  Equation 

Expression  for 
a  System 

Angul ar 
Momentum 

Moment  of  Momentum 

d 
dt^ 

Second  Law 
of 
Thermodynami  cs 

Entropy 

d 

dt 

J 

,s  dV> 
syst 

r 

a  dA 

sys  tT 

Looking  at  the  two  previous  tables,  a  clear  pattern  is  seen 
In  the  form  of  the  expr:'essi  ons  for  the  FUPLAs.   We  may  gener- 
alize this  form  as    r 


d_ 
dt 


syst 


X  dV  =  Y 


(18) 


X  represents  an  extensive  property  that  is  "conserved"  within 
the  system  provi  ded  that  the  right-hand  side  of  Eq.  18  is 
zero  (Y=0). 

Y  is  a  quantity  that  represents  how  the  system  interacts 
with  the  surrounding.   Y's  act  on  the  system  to  change  the  X's. 

Note  that  there  has  been  no  difference  in  our  description 
for  a  system  composed  of  a  differential  quantity  of 


33 


fluid  or,  say,  a  solid  such  as  a  cannon  ball.   We  find  the 
system  description,  however,  inconvenient  in  a  flov/  situation 
because  it  is  difficult  to  visualize  it. 

Suppose  you  can  identify  a  fixed  quantity  of  air  flowing 
into  a  fan  jet  aircraft  engine.   It  would  be  very   difficult  to 
follow  that  identified  mass  of  air  as  it  travels  through  the 
engine  and  experiences  changes  in  properties.   See  Figure  5 
below. 


inlet  I 


compressor 

pOlpiQBnsnBM 


iisd^ 


exhaust 


Figure  5.   A  jet  engine 


It  would  be  much  easier  to  choose  a  convenient  control 
volume,  fixed  with  respect  to  the  engine  for  example  and  then 
investigate  the  changes  in  the  air  properties  as  it  flows  in 
and  out  of  the  control  volume. 

Exercise: 

What  would  be  a  convenient  control  volume  to  describe  what, 
happens  to  air  as  it  passes  through  a  smoking  pipe? 
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SELF-CHECK    QUIZ    NO.     2 
Select    the    proper    ansv/er 

1 .  A    system    i  s : 

a.  the  quantity  of  matter  inside  the  control  volume. 

b.  a  collection  of  matter  of  fixed  identity. 

c.  a  collection  of  mass  of  fixed  position. 

d.  the  set  of  substances  that  makes  up  a  mixture. 

2.  The  Momentum  Equation  is: 

a.  the  expression  of  Newton's  Third  Law  of  Motion. 

b.  the  mathematical  expression  of  a  Law  based  on 
observation  of  nature. 

c.  other  name  given  to  the  Moment  of  Momentum  Theorem 

d.  the  expression  of  the  Law  of  Conservation  of 
Angular  Momentum. 

3.  The  First  Law  of  Thermodynamics: 

a.  cannot  be  applied  to  a  flow  situation. 

b.  has  a  mathematical  expression  usually  called  the 
Entropy  Equation. 

c.  is  applicable  to  the  interchange  of  heat  and  work 
between  the  system  and  the  control  volume. 

d.  has  a  mathematical  expression  called  the  Energy 
Equati  on  . 

4.  The  expressions  for  the  Fundamental  Physical  Laws  for 
a  system: 

a.  show  a  clear  pattern  that  makes  possible  to 
generalize  them  using  a  simple  notation. 

b.  are  difficult  to  apply  in  certain  flow  situations. 

c.  involves  vector  and/or  scalar  quantities. 

d.  all  the  above  are  correct. 


(P)l7  '(P)e  *(q)2  *(q)L   :sJ9Msuv 
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5.   In  the  generalized  expression  for  the  Fundamental 
Physical  Laws  for  the  system  the  element  of  volume 
dV  1s  related  to: 

a .  the  Con  trol  Vol ume . 

b.  a  volume  inside  which  properties  change  from 
point  to  poi  n t . 

c.  the  space  occupied  by  the  system  which  is 
stationary. 

d.  none  of  the  above  is  correct. 


(p)g   :jaMSUv 
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Ill  .   DEVELOPMENT  OF  EXPRESSIONS  FOR  THE  FUMDAME:1TAL  PHYSICAL 
LAWS  USING  THE  CONTROL  VOLUME  APPROACH 

Read  continuously  up  to  the  end  of  Section  B,  but  only  if 
you  are  acquainted  with  the  material  up  to  here.   Otherwise, 
do  not  proceed  until  understanding  the  system  approach.   You 
can  expect  the  reading  of  this  unit  to  take  about  twice  as 
long  as  each  of  the  previous  units. 


A.   THE  DEVELOPMENT  IN  GENERALIZED  NOTATION 

Instead  of  developing  separate  expressions  for  every 
Fundamental  Physical  Law,  all  of  them  can  be  done  at  once, 
using  the  generalized  notation 


r 


d_ 
dt 


X  dV 


(18) 


^ 


syst 


Remember  that  in  the  above  expression 
1 .   Is  valid  for  a 


only. 


2. 


XdV    is    an   extensive    property--a   dependent   quantity 

that    is    conserved   the    system   when    Y    is    zero 

(except    for    the    entropy). 

Y    is    a    quantity   that    affects    the    state    of   the    system. 

The    differential    dV    refers    to    the    volume    filled    by    the 

DO    NOT    CONFUSE    WITH    THE 

RIGOROUSLY    AND    COMPLETELY    DIFFERENT    CONTROL    VOLUME. 

We   will    use    a    control    volume    (CV)    fixed   with    respect    to 
.       Fluid     can     flow   through    this 


3. 

4. 


control    volume.       A   typical    CV    is    depicted    in    Fig.    6. 


Luat^sy^s    8:^euLpjoo3    [en^^UL 
ma:;sX's    jo    ssem    [oj:^uoo 

V  \/  V 

ssem 
ssem    [oj:^uod   jo   m9:;sX's 


'Z 


SjI8MSU\/ 
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cv 

Figure  6.   A  Nozzle  as  an  Example  of  Control  Volume 

I-t  is  possible  to  proceed  with  our  development  using 
the  above  shown  control  volume,  but  a  more  general  form  or 
shape  may  be  preferred,  since  we  are  trying  to  be  general 
in  the  notation.   In  solving  particular  problems  using  the 
control  volume  approach  a  wise  choice  of  a  convenient 
control  volume,  suitable  to  particular  situations,  has  to 
be  made. 
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Let    us    consider    the    following    general    flow: 

At    time    "t"    the    control    volume    is    entirely    occupied    by 
the    system  _(V^^^^    =    V^,^). 


^y^ 


'j^-^  ■y■ 


;^^/ 


^^"  ^^"'y . 


J 


-5^' 


0i^0:^-^':^-':-' 


"\ 


.V^^'.-v'.C-- 


Figure  7.   Snapshot  at  time  "t" 
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At  time  "t  +  At"  we  can  distinguish  three  different 
definite  regions: 


"  ■  '"-'         ■•^!^%^'\:i'<X'7^ 


.-,^---  -"•■'  ■  ;  '  .   '  '-:•;-' .'.i<!S'.;"A'i'.^  ~'"~ '';'---  -Ji?*-- 


.Jar-;<-"' 


':;^^--.=^ 


-^^m^^myim{^e^^^:SM^l^J^'-;;- 


•  V'i^-  -m'^-  -<-^ 


'i^'^:-: 


■•::5^ 


/^; 


:-:^:M'r:''  :S^..i:i:;£:;^:^^'^^^^^^ 


Figure  8.   Snapshot  at  time  "t  +  At' 


Note  that  the  system  may  have  changed  shape  at  time  t  +  At. 
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Region  I:     common  to  the  control  volume  and  the  system 

Region  II:    occupied  by  the  portion  of  the  system  that  left 
the  control  volume  during  the  interval  At 

Region  III:   portion  of  control  volume  unoccupied  by  matter 
belonging  to  the  system  because  it  changed 
position  during  the  interval  At 

From  calculus,  the  time  rate  of  change  of  the  extensive 
property  X  in  the  system  is 


d_ 

dt 


/^ 


syst'^  ^^  =  ^^■^t->0  At 


syst 
@t+At 


X  dV 


X  dV 


(19) 


Remember  that  in  this  case  the  differential  dV  refers  to 
the  volume  occupied  by  the  system  not  the  control  volume 

Now  at  time  t  +  At 

syst  =   region  II  +  region  I 


or 


and  at  time  t 


syst 


syst 


region  II  +  CV  -  region  III 


CV 


(20) 
(21) 
(22) 


Using    these    last    volumetric    considerations 


r 


d_ 
dt 


syst 


X    dV=lim 


At^O 


r 


RII 
(3t  +  At 


XdV  + 


r 


CV 


XdV- 


0t+At         J(at+At 


RIII 


XdV 


CV 

@t 


XdV 


(23) 


Rearranging 

r 


d_ 

dt 


syst 


X    dV=lim 


At^O 


r 


X    dV- 


r 


CV 


X    dV  + 


r 


RII 
J(9t  +  At 


X    dV- 


RIII 

0t+At 


X    dV 


(24) 
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Let    us    v;ork    with    the    right-hand    side    of   Eq.    24    taking    two 
terms    at    a    t i  me . 

The    first    two    terms    represent    changes    with    respect    to    time 
of    the    extensive    property    X    inside    the    control    volume. 


Thus    using    the    definition    of    l_imit: 


lim      1 
At^O    At 


r 


X    dV    -  X    dV      =      -^ 


X    dV 


(25) 


Note  that  the  right-hand  side  of  Eq.  25  has  a  partial 
derivative,  meaning  that  there  is  a  change  with  respect  to 
time  of  a  quantity  that  depends  on  several  independent  vari 
ables.   In  the  control  volume  X  =  X(x,y,z,t). 

We  will  take  care  now  of  the  last  two  terms  of  Eq.  24 


lim  1_ 
At^O  At 


r 


RII 

0  t  +  At 


X  dV  - 


J 


RIII 
Lo  t+At 


X  dV 


These  terms,  associated  with  regions  II  and  III,  represent 
the  amount  of  the  extensive  property  X  inside  these  regions 
at  time  t+At. 

The  amount  of  the  extensive  property  X  inside  II  at  time 
t+At  can  be  related  to  the  amount  of  X  that  can  enter  during 
At  to  this  region,  crossing  that  part  of  the  control  surface 
represented  by  segment  3-4  in  Fig.  9a. 

Similarly,  the  amount  of  property  X  inside  III  at  time 
t+At  can  be  related  to  the  amount  of  X  that  can  enter  to  this 
region  during  the  interval  At  crossing  that  part  of  the  control 
surface  represented  by  segment  1-2  in  Fig.  9a. 

Later  we  will  find  what  this  relation  is,  in  order  to 
transform  the  two  last  terms  in  Eq.  24. 

Referring  now  to  Fig.  9b:   Identify  the  system  drawing  a 
segmented  line  as  in  Fig.  9b.   Also  identify  the  segments  1-2 
and  3-4  of  Fig.  9a.  in  the  blank  circles  in  Fig.  9b.   (Answer 
in  Fig.  B2) 


42 


0 


© 


Control    Vol'jme 
Sys  tern 


Figure   9a.      One    dimensional    version    of   Figure   9b 


Figure   9b.      Two    dimensional    projection 
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Let  us  find  the  mather^ati  cal  expression  of  the  aniounts  of 
property  X  flowing  in  and  out  of  the  control  volume  across  the 
control  surface.   If  dA  is  an  increment  in  the  control  surface 
represented  by  segments  1-2  and  3-4  in  Fig.  9,  then 


V 


n    is    the    component    of    v    perpendicular    to    dA, 
(v    •    n)dA    is    the    volumetric    flow    rate    and 
X(v    •    n)dA    is    the    flow    of   property    X 

Note    that    X(v    •    n)dAAt    =    amount    of    X    that    crosses    dA    in    time    At 

Since    dV    =    dA(v    •    n)At 

Replacing  this  last  expression  for  dV  in  the  last  term  of 
the  right-hand  side  of  Eq.  24 


11m  ]_ 
t^O  At 


r 


r 


RII 
(3t  +  At 


X  dV 


lim  ^_ 

t^O  At 


RII 


X(v  •  n)dA  At 


J  @t+At 

The  integral  in  the  right-hand  side  of  the  above  equation 

represents  the  extensive  property  X  that  flowed  into  region  II 

during  the  interval  At  and  it  was  said  before  that  this  flow  is 

across  that  part  of  the  control  surface  represented  by  3-4  in 

Fig.  9.   Then  defining  A   .,  the  area  corresponding  to  this 
3  ^   out  r      :? 


segment,  one  mav  write 


lim  1_ 
t-^0  At 


r 


y 


RII 

@  t  +  At 

-y 


X(v  •  n)  dAAt 


X(v  •  n)  dA 


(26) 


J   ^out 


Now  if  a  vector  n'  is  defined  as  inward  normal  (inward  and 
outward  refer  to  the  control  surface)  the  quantity 

A.  X(v  •  n')dAAt  will  represent  the  amount  of  property  X 
that  has  flowed  into  region  III  during  At,  so 


lim  1_ 
At^O  At 


r 


r 


J 


RIII 
@t+At 


X  dV  = 


X(v  •  n' )dA 


(27) 


1  n 


and  from  Eqs.  26and  27 

/? 

1  i  m      1/  y     .w 

At->0    Atl  Rir    °^ 

y@t  +  At 


r 


dV 


X(v.n)dA 


J 


out 


n')dA 


(28) 
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but  n    =  -n'  (by  construction) 


lim  1_ 
At^O  At 


RII 
(9t+At 


X  dV  - 


r 

RIII 
^(at+At 


X  dV 


X(v  •  n)dA 


out 

r 


(29) 


J 


A 


X(v  •  n)dA 


in 


-  The  combination  of  the  two  terms  on  the  right-hand  side 
of  Eqs.  28  and  29  represents  the  flow  across  the  boundaries 
of  the  control  volume  that  occupies  regions  III  and  I  at 
time  t  +  At.   We  can  rewrite  Eq.  29  as 


lim  ]_ 
At^O  At 


r 


RII 
(at+At 


X  dV  - 


J 


RIII 
@t+At 


X  dV 


X(v  •  n)dA 


(30) 


Because  the  integration  over  Ain  and  Aout  is  the  same  as  the 
integral  over  the  whole  control  surface  and  nds  e  dA  we  can 
substitute  Eq.  30  into  Eq.  25 


(31) 


r 

r 

f 

d 
dt 

syst 

dV    = 

'6 
9t 

cv^ 

dV 

^( 

^csxt^  • 

dA) 

0 

J 

J 

) 

which  can  also  be  written  as 


r 


d_ 
dt 


r 


^X  dV  = 
syst       ^ 


r 


cv 


II  dV  ^_pesX(v  .  d^) 


(32) 


(Through    the    use    of   Leibnitz    rule    which    is    applied   to   the 
control    volume    in    space.) 

The    first    term   of   the    right-hand    side    may    be    considered 
as    the    rate    of   change    of    the    quantity    of    X    stored    in    the 
control    volume    while    the    other    term    as    the    net    flux    of   X    out 
of   the    control    surface.     '"- 
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Equations  31  and  32  are  very  important  because  they  re- 
late  the  time  rate  of  change  of  property  X  within  a  system 
or  control  mass  with  changes  happening  in  the  control  volume. 

Let  us  examine  the  first  term  of  the  right-hand  side  of 
Eq.  31.   It  represents  the  time  rate  of  change  of  the  exten- 
sive property  X  within  the  control  volume.   It  is  a  triple 
integral  or  a  volume  integral,  so  when  we  write 


8_ 

dt  J 


cv 


X  dV   we  mean 


X  dV 


This  term  is  often  called  "local  time  rate  of  change"  or  storage 

term. 

The  second  term  of  the  right-hand  side  of  Eq.  31  represents 
the  net  efflux  (outflow)  of  the  extensive  property  X  across  the 
boundaries  of  the  control  volume,  also  called  control  surface. 
The  vector  dA  is  perpendicular  to  the  control  surface  at  every 
point  of  it  and  is  positive  when  directed  outward.   The  dot 
product  V  times  dA  gives  the  projection  of  v  over  the  normal 
to  the  surface  direction.   Again: 


r 


o 


X(v  •  dA) 


we  mean 


rr 


X(v  •  dA) 


This  term  is  often  called  "the  convective  term." 

Equations  31  and  32  are  transformation  equations  that  enr 
able  us  to  obtain  an  expression  for  the  FUPLAs  for  a  control 
volume  from  the  corresponding  equation  for  a  system. 

From  Eqs.  31  and  18  we  obtain 


fa'   dv  ^(f) 


9t    CV 


CS 


X  (v  •  dA)  =  Y 


(33) 


This  is  the  generalized  equation  for  the  FUPLAs  in 
control  volume  form. 
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B.   PARTICULARIZATION  OF  THE  GENERAL  EXPRESSIOl 


Now  that  Eq.  33  gives  a  general  expression  for  the 
FUPLAs  in  terms  of  X  and  Y,  it  may  be  applied  to  specific 
cases : 

1  .   Continuity  Equation 

X   =   p 

Y    =    0   /-  r 


or 


or 


at  J 


cv 


pdV 


.o_ 


^ 


cs 


p  V  •  dA 


2 .   Momentum  Equation 

Note  that  here  both  X  and  Y  must  be  vectors. 
t      -      pv 


net  ^ 

It  IcvP^''^  '   ^ 


> 


cs 


pv(v  •  dA)   =   "? 


net 


19  (pv) 
]  9t 


dV 


A 


J 


pv(v  •  dA)  =  t  ^. 
^  net 


CS 


3. 


Energy  Equation 
X   =   pe 
Y   =   Q  -  P 

r  r 

(,y(pe)dV  +^ 


3_ 

at 


cs 


pe(v 


dA) 


(34) 


(35) 


(36) 


=  Q  -  IP     (37) 
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or 


(38) 


In  supplement  B,  Eq.  38  is  transformed  into  a  mo 
convenient  form  by  introducing  the  enthalpy.   This  form 
shown  below  (read  supplement  B  to  see  how  this  was  done) 


re 
i  s 


r 


Q  -y 


r 


''   'Jcv  ft  (p^)  ''  'j 


CS'^'^    2g    g  ^'^^^ 
^c   ^c 


dA) 
(39) 


Where  ]Ppw  is  the  net  useful  (shaft)  work. 


CV 
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SELF-CHECK    QUIZ    NO.     3 


Select    the    proper    ansv;er 

1.       In    the    generalized    notation    expression    for   the    FUPLAs 


r 


d_ 
dt 


XdV 


The 

a. 
b. 
c. 
d. 

At 
a . 

b. 
c . 
d. 

Whi 
on 


a. 


sysf 

quanti  ty    X : 

must    be    a    scalar. 
is    the    unknown^of   the    problem, 
represents    an    extensive    property, 
none    of   the    above. 

time    "t"    the    control    volume    and   the    system: 

are   held   stationary   with    respect    to    an    inertial 

frame    of    reference. 

occupies    momentarily    the    same    space. 

are    partially    contained    one    inside    the   other. 

cannot   exist   at    the    same    time. 

ch    of   the    following    expressions    is    correct    based 
calculus    (noting    that    in    the    CV,    X    =    X(x,y,z,t)) 


b. 


c . 


d. 


r 


lim      ]_ 

At->0    At 


lim      1_ 

At^    At 


lim      1_ 

At->0    At 


lim      1_ 
At^O    At 


r 


syst^dV   =    ^  J 


CV 


XdV 


=    ^  XdV 


=    ~  XdV 


•(3)C    '(q)2    '(^)L       :sJ8MSUV 
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4.   The  following  equation  relates  the  expression  for 
FUPLAs  for  a  system  with  an  expression  valid  for 
the  con t ro  1  vo 1 ume 


dt  Jsyst^^^    9t  J 


,XdV  + 


cv^"'  jr CS 


9 


X(v  •  dA) 


This  equation  is  directly  derived  from  observa- 

ti  on  of  nature  . 

The    first    term   of    the    right-hand   side    is    really 

a    volume    integral    and    is    known    as    the    "convective 

term. " 

The  second  term  of  the  right-hand  side  of  the 

equation  is  a  surface  integral  sometimes  called 

the  "convective  term." 

Can  be  applied  only  to  unsteady  flows. 


•  (D)|7    :J8MSUV 
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IV.   EXTENSION  QF^THE  CONTROL  VOLUME  APPROACH  TO  OTHER 

FUNDAMENTAL  PHYSICAL  LAWS 


The  Law  of  Angular  Momentum  and  the  Second  Law  of  Thermo 
dynamics  for  a  system  were  stated  without  any  further 
considerations  in  Table  III. 

Now  we  will  try  to  obtain  expressions  of  these  laws 
applicable  to  a  control  volume. 


A.   ANGULAR  MOMENTUM 

From  the  development  of  the  momentum  theorem  for  a  system 
in  Table  II  or  Eq.  12 


/^ 


net 


d 

dt 


syst 


V  p  dV 


(40) 


Performing    the    cross    product    of    both    sides    of    the    above 
equations    with    a    vector    r   whose    origin    is    fixed   to   the    same 
inertial    frame    of   reference    relative    to   which    v    is    measured, 
we    can   write 


net 


d_ 

dt 


/^ 


syst 


r>tp   v 


dV 


(41) 


The    quantity    in    the    left-hand    side    of   the    above    equation 

is    the   moment   of   the    force    F      .    with    respect    to   the    origin 

net         ^        ^     ^ 

of  r.   We  have  assigned  to  the  moment  the  symbol  M^^. . 

^  ■^  net 

The  integrand  on  the  right-hand  side  of  Eq.  41  is  the 

cross  product  of  the  moment  arm  r  and  the  momentum  and  defines 

a  quantity  called  moment  of  momentum  or  angular  momentum. 

Using  the  generalized  notation  introduced  previously: 


p  (  r  X  v) 


t 


^ 


net  .. 


(42) 
(43) 
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Figure  10.   Coordinate  System 


Using  the  expression  for  the  Fundamental  Physical  Laws 
develop^^d  for  a  control  volume  in  Eq.  29 

r 


r 


r 


d_ 
dt 


So 


J 


syst^^"  X  v)dV  =  ^ 


CV 


p( r  X  v) 


v)dV+^'^. 


J 


CS 


p(  r  X  v) ( V  dA) 


(44) 


This  is  the  so-called  Moment  of  Momentum  equation. 

B.   SECOND  LAW  OF  THERMODYNAMICS 

The  Entropy  Equation  for  a  system  as  the  mathematical 
expression  of  the  Second  Law  of  Thermodynamics  was  stated  in 
Table  III. 
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inequality  valid  for  an  irreversible  process 


4  Ids 


(45) 


equality  Vclid  for  a  reversible  process 


Where  6Q  is  a  small  quantity  of  heat  transferred  during 
an  elementary  part  of  the  cycle,  T  is  the  absolute  tempera- 
ture at  that  point  of  the  boundary,  S  the  entropy,  a 
property  and  s  the  specific  entropy  (entropy  per  unit  mass) 

r 


syst 


spdV 


(46) 


Also 


6Q 


r 


area 


qdt  dA 


(47) 


where  q  is  the  heat  flux  (amount  of  heat  flowing)  per  unit 
time  per  unit  area 


dA  ^dt^ 


(48) 


Dividing  both  sides  of  Eq.  47  by  T 


r 


6Q 
T 


%   dt  dA 
area  T 


(49) 


Substituting  Eq.  45  into  Eq.  49 


but  from  Eq.  46 


dS  >   I     3.  dt  dA 
—  Jarea  T 


r 


dl  > 
dt  - 


J 


%   dA 
area  T 


r 


dS^  ^   d_ 

dt     dt  •^syst 

and  finally  from  Eqs.  51  and  52 


d_ 
dt 


sp  dV 


/ 


'J 


syst 


spdV  > 


^ 


area  3-  dA 


(50) 
(51) 


(52) 


(53) 
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C.   SELF-CHECK  QUIZ  NO.  4 

1.  Can  you  associate  Eq.  53  to  the  generalized 
notation  that  we  had  already  developed? 

2.  What  is  the  difference,  besides  the  type  of 
quantities  involved,  between  Eq.  53  and  the 
expressions  for  the  other  FUPLAs  for  a  system 
in  Table  II? 

3.  Can  you  develop  an  expression  for  the  Second  Law 
of  Thermodynamics  valid  for  a  control  volume? 

(See  answers  on  next  page) 
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Answers  to  Self-Check  Quiz  No.  4 


1  .   X  =  PS 

r 


Y  > 


J 


A     T-  dA 
Area  T 


2.  The  sign  >.  depends  on  what  type  of  process  is  taking 
place  for  the  heat  being  transferred  in  or  out  of 
the  system.   The  equal  sign  holds  if  the  process  is 
reversible.   The  inequality  sign  holds  if  the 
process  is  irreversible. 

3.  Applying  X  and  Y  to  Eq.  29 


T 


8_ 

at 


r 


J 


^^PSUV    Jj^QS 


r 


dV  +(!)pcPs(^^  •  dA) 


J 


area  T 


dA 
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V.   APPLICATION  TO  PARTICULAR  PROBLEMS 


A.   CONTINUITY  EQUATION 
From  Eq .  34 . 

r  r 


J 


cv 


iP 
9t 


dV 


.  C) 


JCS 


p(\r  •  dA) 


Applying  this  equation  to  the  following  case: 

One-dimensional,  steady,  incompressible  flow  as  shown  in 
Fig.  11. 


Figure  11.   Container  with  unequal  inlet  and  exit  areas 


One-dimensional  flow  means  that  the  velocity  is  a  function 

of  coordinate.   There  is  no  change  of 

velocity  with  respect  to  other  coordinates. 

Steady  flow  means  that  the  velocity  is  not  a  function  of 


[B  Lt^eds  8U0  R\_  uo 
-.sjawsuy 
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Incompressible  flow  means  that  the 
change . 


does  not 


Considering  the  above,  we  can  immediately  make  the  first 
term  of  Eq.  34  equal  to  zero. 


Now  Eq.  34  becomes 

r 


J 


cs 


p  V  .  dA   =   0 


(53) 


p  is  a  constant  so  we  can  take  it  out  of  the  integral  in 
Eq.  54. 

Also  V  is  a  function  of  x  only  so 


V 


V  1 
X 


And 


But 


dA   =   dA  t  +  dA   J  +  dA   k 
X       y  ^     z 


V  •  dA  =  V  dA   +  0  dA   +  0  dA 
XX        y        z 


So  Eq .  54  becomes 

r 


V  dA 

X    X 


CV 


(54) 


It  is  only  necessary  to  evaluate  the  integral  Eq.  54  to 
stations  1  and  2  of  the  control  volume  since  there  is  no  flow 
through  the  wal Is 


(P  \^x\    -    (P  \^x\ 


0 


Eliminating  the  subscript  x:   v   =v-,    and   v   =v 

X  ^      I  o 


v^  A^   =   V2  A^ 


(55) 


In  the  evaluation  of  the  integral  e'quation  care  must  be 
taken  with  the  signs  remembering  that  dA   can  be  positive  or 

A 


/Ci^LSuap      :jaMSU\/ 
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negative,  depending  on  the  

(i.e.,  at  the  inlet  v  •  dA  is  a 


u  a  n  t  i  t  y ) 


In  the  following  example  we  will  show  the  system  and  con- 
trol volume  approach  to  an  unsteady  problem. 

B.   EXAMPLE  OF  APPLICATION  OF  THE  SYSTEM  AND  CONTROL  VOLUME 
APPROACH  TO  THE  SOLUTION  OF  AN  UNSTEADY  PROBLEM 

A  partially  evacuated  bottle  initially  containing  a  mass 
m,  is  to  be  filled  with  air  from  the  atmosphere  until  the  mass 
reaches  a  value  of  m^ .   Find  the  work  to  be  done  to  achieve 
this  and  if  there  is  no  heat  interaction,  find  the  change  in 
internal  energy. 

1  .   System  Solution 

A  system  or  control  mass  is  . 


We  have  to  choose  a  system  consistent  with  the  above 
definition.   The  mass  during  the  process  must  remain  unchanged. 
For  this  purpose  we  define  the  system  boundaries  so  as  to 
enclose  the  final  mass  mp  at  all  times.   (See  Figures  12  and 
13). 

We  are  consistent  with  the  definition  because  this 
system  contains  a  fixed  quantity  of  mass  (m^)  although  the 
volume  that  contains  it  changes. 

The  boundaries  of  the  system  at  time  t  contain  mass  m-, 
inside  the  bottle  and  mass  m^  -  m,  outside  the  bottle.   The  mass 
of  fixed  identity,  i.e.,  the  system  is  


(write  an  algebraic  expression  in  terms  of  m,  and  m^) 


'■01     +     ('■LU     -     ^LLl)      =     ^lU 

A:^L:^U9pL    paxL^    ^o    sseiu   j.o    A:iL:^uenb    g 

9AL:^e6au 
BauB    9L|:;    6  u  L:^u8sajdaj    jo;d8a    8lj:^   j.o    uoLttoauip    aij:; 

:  SJ8MSU\/ 
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^  m^ -m. 


'•^ .^ 


system  - 


insulated 
bottle 


Figure  12.   System  solution.   State  1 


insulated 
bottle 


system 


Figure  13.   System  solution.   State  2 
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At  time  t  +  At  the  boundaries  of  the  system  have 
changed  forming  a  new  volume,  but  being  consistent  with  the 
definition  of  system,  the  mass  contained  in  it  is  still 


Using  the  First  Law  of  Thermodynamics  with  the  sign 
convention  of  Fig.  4  we  may  write 


1^2 


1^^2 


=      E^      -      E^  '  (56) 

The    first    term   on    the    left-hand    side    of   Eq.    56    may    be 

eliminated  because 

If  the  only  changes  in  energy  are  changes  in  internal 
energy,  we  can  apply  the  above  mentioned  expression  of  the 
First  Law  of  Thermodynamics  for  a  system,  reducing  it  to 


l"2 


"2  -  "l 


The  change  in  internal  energy  is  the  internal  energy 
at  t+  At  minus  the  internal  energy  at  t 


AU   =  m^u^    -    [m^u^  +  (m^  -  m^ )  u^^J 


(57) 


Where  the  subscript  "atm"  refers  to  the  properties 
of  the  air  at  atmospheric  conditions  outside  the  bottle.   The 
atmosphere  acts  as  a  reservoir  (all  properties  remain  constant) 

The  work  -.W^  is  the  work  that  has  to  be  done  on  that 
part  of  the  system  outside  of  the  bottle  at  time  "t"  to 
introduce  it  to  the  bottle. 


l"2 


p  .   1^  ^  [0  -  (m^  -  m,  )  1 
^atm   atm"-     ^2     r  -^ 


(58) 


and  Eq.  57  becomes 


'^atm   atm  ^2   1 


m 


2^2-^^^l  -  ('^2-'^l^  ^tm^ 


(59) 


2 .   Control  Volume  Solution 

If  we  now  solve  the  same  problem  using  the  control 
volume  technique,  we  have  to  choose  the  control  volume 
conveniently. 


pa:;e[nsuL  sl  3[:;:^oq  b[\i 

Z 


LU    :SJ9MSUV 
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insulated. 

b  0 1 1 U^ 


Patm  "atm 


CV 


Atmosphere 


Figure  14.   Control  volume  solution 


During  the  process  of  filling  the  bottle  the  atmos- 
phere is  a  , 

^tm^  


Its  properties  (i.e.,  P^^^,  u^^^  and 


Recalling  Eq.  39,  the  general  expression  for  the 
Energy  Equation  for  a  control  volume. 

r  r 

^c      -"c 


Q  -  P 


CV 


_(ep)  dV  +  Y  (h  +  i^  +  z  S_)p(^.dA) 

cv^       Jcs    ^9^     g, 


In    the    problem    under    consideration    there    is    no    heat 
transfer   and   we    can    neglect    both  p^w    and   the    potential    energy 
changes.      Thus    Eq.    39    reduces    to: 


fv-(ep)    dV    +Y       ^^^   ^Yi-)p(v-dA) 
CV  ->'CS  'C 


0 


(60) 


We  are  solving  a  non-steady  flow  problem.   Integrating 
Eq.  60  with  respect  to  time,  t^  being  the  time  at  which  the 
bottle  is  filled  with  mass  m„ ,  we  have 


paBupqoun  uLeiuaj 

jLOAjesaj   :sj9Msu\/ 


61 


J 


|-  (ep)dV  +(^ 

cv  Jcs 


0 


[h     +  ir^  )  P  ( V  •  dA) 


2g, 


dt   =   0    (31) 


Working  the  integral  Eq.  G 1  by  terms,  we  have  for  the 
first  term 


A 


^ 


0 


r 


y 


cv 


.   -9T  (ep)  dV 


dt 


/^2 


1 


dE 


(62) 


And  if  all  changes  in  energy  inside  the  Control  volume  are 
changes  in  internal  energy,  Eq.  62  reduces  to: 


^t. 


0 


r 


L- 


J 


cv 


3t 


(ep)  dV 


dt 


H  -   ^1 


(63) 


Where  the  subscripts  2  and  1  mean  final  and  initial  states 
respectively.   Now  let  us  define 


h.   E  h  +  ;L^ 

t       2g, 


(64) 


the  total  or  stagnation  enthalpy.   Replacing  Eq.  64  in  the 
second  term  of  Eq.  61  we  obtain: 


0 


r 


-\ 


J 


CS  (h  '    2^)p(v-dA) 


dt 


rt. 


J   0 


r 


J 


CS 


h^  p(v-d^) 


dt 


(65) 


dm 


But  p(vdA)  =  -^   .       For  details  on  this  see  Fig.  9  and  the 


development  of  Eq.  26.   With  this,  Eq.  65  reduces  to: 


-</ 


0   L 


Z' 


y 


CS 


(h  +  ^)p(v.dA) 


dt 


/•t. 


>  0 


/.  •  dmx    u  dmx 


1 


t  dt^      dt 

out      in 


dt 
(66) 
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There  is  no  outward  flov;  of  mass.   The  total  inflow  of 
mass  during  the  interval  from  time  0  to  t^  is  (rip  -  m-,  )  . 
During  this  interval,  the  total  or  stagnation  enthalpy  of  the 
air  entering  the  control  surface  is  equal  to  the  constant 
atmospheric  enthalpy  (h  o^  ).   With  all  these  considerations, 
Eq .  66  reduces  to : 

r 


rt. 


0 


(h  +  ^)p(^.dA) 
CS      ^^c 


dt   =   -  (m^  -  m^)  h^^^    (67) 


Substituting  Eqs.  63  and  67  in  Eq.  61  we  obtain: 


And  if 


^2  -  ^1   =   ^tm('^2  -  ^1^ 


Up    =    Up  Hp 


atm 


atm   atm 


""•^at 


m 


1  '"1 


Patm  'l^atm^-'Z  "  "'1 


(m^  -  m^)   =  m^u^    -    [m^u^  -  (m^  -  m^)   u^^^^ 


(68) 


Equation  68  is  exactly  the  same  as  Eq.  59  obtained  using 
the  system  approach  even  though  a  few  additional  steps  were 
needed. 
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VI.   SUPPLEMENT  A:   BODY  AND  SURFACE  FORCES 


In  discussing  the  Momentum  Equation  we  had  in  Eq.  11 

net     /      volume  ,  /      surface 

And  at  that  time  the  distinction  between  a  body  and  a  surface 
force  was  not  explained  in  detail. 

The  body  or  volume  forces  are  proportional  to  either  the 

mass  or  volume  of  the  body  and  comprise  those  forces  involving 

action  at  a  distance.   These  forces  are  a  consequence  of  long 

range  force  fields,  and  include  such  examples  as  the  force  of  . 

gravitational  attraction,  magnetic  forces  and  electrical 

(coulomb)  forces.   For  accelerating  coordinate  systems,  "inertia 

forces"  such  as  centrifugal  and  coriolis  forces  would  also  be 

included.   But  in  this  case  these  last  two  forces  must  not  be 

consicdered  because  our  development  was  done  for  an  inertial 

2 
coordinate  system. 

Surface  forces  are  those  forces  which  are  exerted  at  the 
control  surface  by  the  material  outside  the  control  volume  on 
the  material  inside  the  control  volume.   These  forces  are  a 
consequence  of  short  range  force  fields.   Such  forces  are 
exerted  in  the  form  of  surface  stresses.   We  can  distinguish 
two  types  of  surface  forces:   (1)  those  arising  from  normal 
stresses,  or  pressures,  acting  on  the  control  surface, and 
(2)  those  arising  from  shear  stresses,  or  viscous  stresses, 
acting  on  the  control  surface. 


See  footnote  on  page  22 
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VII  .   SUPPLEMENT  3 


INTERMEDIATE  STEPS  IN  THE  DERIVATION 


OF  THE  ENERGY  EQUATION 

Applying  the  general  expression  for  the  FUPLAs.  Eq.  31, 
to  the  Energy  Equation 


r 


cv 


9(pe) 
8t 


dV  +()   pe(v-dA) 
J   CS 


Q 


iP 


(38) 


Also  stated  without  proof  was  that  Eq.  38  can  be  transformed 
into  Eq.  39.   To  do  this,  what  F  is  must  be  considered.   In 
Eq .  16,  which  is  the  most  elementary  expression  for  the  First 
Law  of  Thermodynamics,  was  stated  that  P  was  the  rate  of  work 
done  by  the  system  to  the  surroundings. 

In  Supplement  A  the  nature  of  the  forces  that  may  be  present 
and  dividing  them  into  body  and  surface  forces,  was  examined. 
Body  or  volume  forces  do  not  appear  in  P  since  they  do  not 
represent  work  done  on  the  surroundings.   Rather,  they  appear 
in  the  energy  term  as  a  potential  energy  (for  most  cases  of 
interest).   We  write  e  as: 


u  + 


v|2 


gz  + 


Hence,  we  restrict  ourselves  to  the  contribution  of  the  s  urf ace 

forces  to  IP.   These  comprise  normal  and  shear  stresses  and  the 

objective  here  is  to  single  out  the  contribution  of  the  normal 
stresses  which  are  assumed  to  be  equal  to  the  static  pressure 
acting  on  the  fluid. 

Consider  an  element  of  fluid  of  mass  pdV  entering  the 
control  volume. 
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dA 


F.igure    Bl  .       Diagram    for    calculating    flow    work 


work 


rate  of  work 


force  X  displacement 


(pdA)dX 


1 


(pdA    dX) 


1 


(pdV) 


At    ^^""    ^"'  At 

Hence    dividing   by    the    mass    (pdV)    we    get 


rate  of  work/lbm 


k  (P/P^      =      It    (P*^) 


(Bl) 
(B2) 

(B3) 


The    term      p/p    is    called    flow   work    or    flow    energy    since    it 
represents    the   work    to    get    the    fluid    into   or   out    of   the    con- 
trol   volume.      The    important    thing    to   note    is    that   we    can 
separate    out    the    contribution    of   the    normal    stresses    to  P   by 
simply    accounting    for   the    flow   work.      Since    the    fluid   must 
always    perform    this    work    upon    entering    or    leaving    the    control 
volume,    this    is    a    very    useful    step. 


T^E    1/p    is    the    specific   volume    (volume    per    unit   mass, 
3 
i.e.,    ft    /Lb    ).       Do    not    confuse    this    symbol  O^with    V    for 

volume    or   v    for   velocity. 
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IP 


pw  +  (rate  of  flow  work) 


(B4) 


In  Eq.  B4,  Ppw  represents  what  remai ns,  namely,  work  done 
by  shear  stresses  and  shaft  work.   Since  friction  is  usually 
handled  indirectly,  Ppw  is  commonly  labeled  as  the  net  useful 
or  shaft  rate  of  work. 

Let  us  see  what  the  above  means  for  a  simple  compressible 
substance.   Suppose  a  fluid  enters  the  control  volume  at 
station  (1)  and  exits  at  station  (2)  and  suppose  further  that 
all  the  properties  are  homogeneous  at  each  station.   Now,  on 
a  per^pound  basis,  the  fluid  must  do  an  amount  of  work  P  i  l^i 
to  enter  the  control  volume,  and  an  amount  of  work  Po  i^o  to 
leave  the  control  volume.   The  work  that  may  be  performed  is 
still 


1^2 


r2 


J  1 


r2 


p  d^   = 


Vdp 
P 


so  that  analogous  to  Eq.  B4, 

w 


Pl^l^ 


dp 


(this  is  an  important  result  which  is  usually  written  down 
without  proof.   Show  thi  s  to  yoursel f  starting  with  a  steady 
flow  form  of  the  First  Law  and  setting  AKE  =  0  and  APE  =  0). 

Now  we  are  ready  to  arrive  at  the  final  form  of  the  First 
Law  namely,  Eq.  39.   Referring  back  to  Fig.  9  we  can  do  this. 
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Figure    B2.      Two-Di men s i onal    Projection    at    time    t+At 


Equation    B2    can    be    written    as 


1^  (p    dAdx)      =      ^  (p    dA    •    V   At) 


(B7) 


Since  for  the  fluid  entering  dA  opposes  v  and  for  the  fluid 
exiting,  they  are  in  the  same  direction,  we  get 


rate  of  flow  work 


lim   ]_ 

At^O  1  At  1 Sout 

1 


p(v'dA) 


At 


Si  n 


P( 


v-d/^)  J 


=  O    p(v-dA) 
J  CS 
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(B8) 
(B9) 


Hence,  from  Eq.  B4 


^    =    ^a  ^ 


cs 


(v-d^) 


(BIO) 


And  the  First  Law  becomes 

r  r 

3_(_eel  dv  +  C)   [u  +  ^  + 
^CV    dt       JcS      P 


yL+  za_]p(^.dA) 


Q  -  IP 


CV 
(BID 


Now  substitute  h  =  u  +  ^  to  get  Eq.  39. 
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VIII.   SUPPLEMENT  C:   EXPRESSIONS  FOR  THE  FUNDAMENTAL  PHYSICAL 

LAWS  (FUPLAs)  IN  A  DIFFERENTIAL  FORM 


The  equations  derived  in  the  main  text  of  the  handout  are 
in  integral  form.   There  are  many  applications  where  it  is  more 
useful  to  begin  with  a  differential  form  of  the  equations. 

An  alternate  for  Eq.  34,  the  general  expression  for  the 
FUPLAs  is 


3t 


dV  +( 


O 


y 


cs 


where 


X(v  •  dA) 


y  dV 


r 


J  cv 


y  dV 


(CD 


provided  Y  is  expressible  on  a  per  unit  volume  basis. 

Reference   14  or  any  other  Vector  Analysis  book  explains 
in  detail  a  very  useful  mathematical  tool,  the  Divergence 
Theorem. 

The  Divergence  Theorem  (also  called  Theorem  of  Gauss) 

->- 
states  that,  for  any  volume  in  a  vector  field  c  ,  the  normal 

component  t * n   integrated  over  the  area  is  equal  to  the 

divergence  of  this  vector  field  integrated  over  the  volume. 

In  mathematical  symbols  we  have: 


t-n  dA 


r 


J  V 


(V-c)dV 


Reference   3  states  that  c  n^ay  be  generalized  to  any 
tensor  field.   Now  let  us  apply  the  Divergence  Theorem  to  the 
second  term  of  Eq .  CI , 


r 


X(v  .  dA) 


V-(X  v)dV 


(C3) 


CV 
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where  it  is  not  difficult  to  recognize  that  we  have  identified 
X    in  Eq.  C2  with  Xv  .   Note  that  this  quantit.v  is  a  vector  if 
X  is  a  scalar  or  a  tensor  if  X  is  a  vector! 


Replacing  Eq.  C3  in  CI 


9t 


X  dV   + 


r 


V-X  V  dV 


cv 


->'cv 


y  dV 


(C4) 


y 


cv 


Therefore,  since  the  control  volume  is  arbitrary  and  common 
to  al 1  three  terms 


91 

at 


+   V-(Xv) 


(C5) 


This  is  the  differential  form  of  the  Fundamental  Physical 
Laws  using  the  generalized  notation  X.   We  will  apply  Eq.  C5 
toentriesinTablell. 


A.   CONTINUITY 


p   and   y 


1^  +  V-pv   =  0 


For  steady  flow 

V-pv   =   0 

4      / 
and  for  incompressible   flow  (p 

Thus 

V'V   =   0 


(C6) 


(C7) 


(C8) 


Some  books  (like  Ref.  4,  p.  44)  goes  straight  from  Eq 
06  to  C8  without  assuming  steady  flow.  'This  can  be  done  if 
the  flow  is  incompressible,  since  then  p  does  not  change 
either  with  time  or  position. 
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OMENTUM 


pv 


and 


y  =  f 


net 


^^net  '^^^  ^"""^  volume) 


--^     +   V.(pvv)   =   f 


net 


(C9) 


Note  that  the  product  vv  is  not  defined  in  vector  analysis. 
It  is  a  dyadic  product  (a  tensor)  and  we  should  be  using  tensor 
notation.   However,  V»(pvv)  represents  a  vector  [Ref.  15,  p. 
200]. 


The  first  term  of  Eq.  C9  can  be  written  as 


9(pv) 

at 


v^  + 


at 


a_v 
at 


and  the  second  term   can  be  written 


(CIO) 


V* ( pvv ) 


v(  V*  pv)  +  (  pv  V)  V 


(Cll) 


Replacing  Eqs .  CIO  and  Cll  in  Eq.  C9 


V  (|^  .   V.pv)   .   (pv.V)v  .  p|^  =   f^^^ 


(C12) 


The  parenthesis  of  the  first  term  in  Eq.  C12  is  zero  from 
continuity  Eq.  C6,  so  Eq.  C12  reduces  to 


av 


->   >  -> 


p[-9Y  +   (v.V)v] 


f 


net 


(C13) 


There  is  a  compact  notation  often  used  in  connection  to  these 
types  of  equations 


2L1  .  iLl  .  (^.v)c  ] 


(C14) 


^Ll  -js  called  the  "substantial,"  "total,"  or  "flow"  deriva- 
tive according  to  different  authors.   Rewriting  Eq.  C13  using 
the  notation  introduced  in  C14 
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d  V     "^ 

which  is  an  equation  of  motion  analogous  t o  m  — -  =  F   . 

dt  net 

the    expression    for    Newton's    second    law    for    a    solid, 
such    as    a    cannon    ball.       A    deeper    insight    is    necessary 

regarding    the    term    on    the    right-hand    side    of    Eq.    C15.      This 

->■ 

term  is  related  to  F   .  whose  meaning  was  explained  in 

net  3       r 

supplement  A.   So  we  can  write 


f 


net 


f   T     +  f    x 
volume     surface 


(C16) 


For  the  moment  we  will  concentrate  our  discussion  on 

surface  since  it  is  not  immediately  apparent  how  to  express 
it  on  a  per  unit  volume  basis.   The  surface  forces  depend  on 
the  rate  at  which  the  fluid  is  strained  by  the  velocity  field 
The  system  of  forces  determine  a  state  of  stress  and  the  re- 
lationship between  stress  and  strain  can  only  be  given 
empirically.   Recalling  again  Supplement  A,  a  distinction 
can  be  made  between  forces  arising  from  normal  stresses  and 
forces  arising  from  shear  stresses. 

f.     =f     i+fu  (C17) 

surface      normal     shear        ^    ' 

The  force  due  to  the  normal  stress  may  be  written: 

f     T   =   -  V  p  (C18) 

normal         ^  ^    ' 

and  the  force  due  to  the  shear  stress,  under  the  Stokes 

hypothesis, 

■f 
shear 

Where  u  is  the  viscosity  of  the  fluid. 

the  fluid  (and  may  be  strongly  dependent  on  temperature) 

Replacing  Eqs.  C19  and  C18  in  C17 


-t  2  ^ 

f_u =   y  V   V 


(C19) 
It  is  a  property  of 


f. 


-  V  p  +  y  V   V 


(C20) 


surface 

We  have  only  shown  between  Eqs.  C16  and  C20  how  the  different 

-> 

kinds  of  stresses  make  up  f    ,.    .  A  complete  derivation 

^   surface  ^ 

of  Eq.  C20  can  be  found  in  Ref.  4.   In  the  derivation  of  Eq. 

C20  it  was  assumed  that  the  fluid  is  isotropic,  namely,  the 

components  of  the  stresses  are  equal  in  all  direction  and 

that  it  is  Newtonian,  i.e.,  a  linear  relation  exists  between 
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the  rate  of  strain  and  the  stress,  the  relation  between  th*^ 
two  being  empirical. 

Replacing  Eq.  C20  in  CI  6 


f.=fT  -Vp+yVv 

net  volume  ^        ^ 

and  replacing  Eq.  C21  in  C15 


(C21) 


(C22) 


Equation  C22  is  a  very    useful  form  of  the  equations  of  motion 
called  the  Navier-Stokes  Equations  (a  set  of  three  equations, 
one  for  each  coordinate).   They  form  the  starting  point  for  a 
course  in  boundary  layer  theory.   Written  as  in  Eq.  C22  it 
has  only  the  restrictions  of  isotropic,  incompressible,  New- 
tonian flows.   The  incompressible  restriction  arises  from  the 
fact  that  we  had  treated  y  as  a  constant  taking  it  out  of  the 
Laplacian  in  the  second  term  of  the  right-hand  side  of  Eq. 
C22.   For  compressible  flow,  where  there  exists  gradients  of 
temperature  that  affect  strongly  the  value  of  p ,  the  above 
mentioned  term  should  be  written  V»yVv.   In  the  derivation  of 
Eq.  C22,  the  compact  form  of  Navier-Stokes  equations,  most  of 
the  effort  was  given  to  explaining  the  significance  of  each 
term  and  its  origin.   For  a  complete  and  rigorous  derivation 
of  these  equations  we  suggest  the  use  of  Ref.  4  in  the  con- 
text of  an  appropriate  course. 

In  Aeronautical  Engineering,  the  Navier-Stokes  equations 
are  useful  in  flows  (not  fluids)  where  the  viscous  effects  are 
important,  i.e.,  flows  where  large  velocity  gradient  exist. 
Such  situations  are  found  in  boundary  layers  where  there  is  a 
change  of  velocity  between  a  value  of  zero  at  the  wall  and  a 
value  different  from  zero  at  the  free  stream.   The  gradient 
creates  important  viscous  forces  taken  i*nto  account  in  the 
last  term  of  Eq.  C22.   In  the  free  stream  outside  the  boundary 
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layer,  although   -(the  viscosity)  exists  as  in  the  boundary 
layer  region  (same  fluid),  the  viscous  forces  do  not  exist 
because  there  is  no  velocity  gradient,  and  the  last  term  in 
Eq.  C22  is  zero. 

In  which  other  case  would  you  do  the  same  ?  

In  these  two  cases  Eq.  C22  reduces  to 


Dt 


f   ,     -  Vp 
vo 1 ume     ^ 


(C23) 


Write    Euler's    Equation    of  Motion    for   one    dimensional    flow 
using    the    definition    of   total    derivative 

(C24) 


For  steady  flow  and  neglecting  the  volume  forces,  Eq. 
C24  reduces  to: 


dv 


1  d£. 
"  p  dx 


dx   -   -  -  ^  (C25) 

Why  does  Eq.  C25  only  contain  total  derivatives? 


Now  integrating  Eq.  C25  we  get 


r 


d£ 

P 


Or  if  the  flow  is 
be  integrated  to  give 

v2 
2 


P 


const  an  t 


constant 


(C26) 


,  the  pressure  term  can 


(C27) 


(:;ue:;suoD    sl    d)    a  [  q  l  ssajdmoou  l 
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This  is  called  Bernoulli's  Equation.   Sometimes  it  has  an 
g 
additional  term   ^g    coming  from  the  body  forces  -.hen  vie 

neglected  in  Eq.  C24  (as  they  pertain  to  a  gravity  field). 

This  equation  looks  like  an  Energy  Equation  but  it  is 
not  complete  as  we  shall  see  in  part  C  of  this  Supplement. 

Exe  rci  se : 

What    does    it    mean   when    the    flow    is    Newtonian? 

On    what    property    does    the    viscosity    primarily    depend? 

It    can    be    shown    that    the    assumption    of   one-dimensional 
flow   can    be    replaced    by   the    restriction    that    the    constant    in 
Eq.    C27   takes    different    values    for   different    streamlines. 
In    other   words,    the    left-hand    side    of    Eq.    C27    is    equal    to    a 
different    constant,    depending    on    which    streamline    of   the    flow 
is    considered. 

Also,    it    can    be    shown    that    Eq.    C27    is    valid    everywhere 
(no   matter    if   we    change    streamline)    if   the    flow    is 
irrotational    (Vxv=0). 


Equation    C27    is    known    as 


Equation.      The 


assumptions    and   constraints    in    its    derivation    are: 


•p9:iD8L68u    S8DJ0J.    aiun[OA      •pLn[j.    ppsLAui 
MO  Lj.    [e  uo  L  suaai  Lp- auo    'Apeat^s    'a  [  q  issajdmoDUi 

L  [  L  n  0  u  J  a  a 
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C.   ENERGY 

In  our  development  of  the  Energy  Equation  for  a  control 
volume,  Eq.  38  we  stated  that  this  form  was  not  yery    con- 
venient and  in  Supplement  B  we  transformed  it  into  Eq.  39. 

The  differential  form  of  the  Energy  Equation,  Eqs.  38 
or  39  can  be  derived  in  a  manner  similar  to  the  Momentum 
Equation,  thus  only  the  results  will  be  given  here.   For  a 
description  of  the  various  terms  and  their  application  see 
Ref.  5. 

The  following  form  is  usually  quoted  [Ref.  5]  for  the 
Energy  Equation: 


p||  +   p(v-V)e   -   pV*v   =   -V»q   +  y$ 
or  in  terms  of  the  enthalpy 

p|y  +  p(v-V)h  +  [||-  +  (v-V)p]  =  -V-q  +  y-D 

where   e 

Q   =   VK   q 


(C29) 


(C30) 


u  (neglecting  changes  in  kinetic 

^-jT   .  ^  and  potential  energy) 

U       Q.dA  ,   ,^ 

^  C31 

Jcs 


and      $   =      dissipation    function    acting   through    the    fluid 
viscosity 

If   neglecting    radiation    heat    transfer   and    assuming    that 
Fourier's    Law   of   Heat    Conduction    applies,    then 


q 


-  k  V  T 


and 


V'q 


V'q 


-  V(k  V  T) 


-  k  V 


for  k  =  constant 


(C32) 
(C33) 
(C34) 


The  solution  to  the  equation  above  represents  the  effect  of 
fluid  motion  on  heat  transfer  and  is  studied  under  the 
title  of  "Convective  Heat  Transfer." 
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Substituting  Eq.  C34  into  Eq.  31 


Pj^+    p(v.V)h  +  [||  +  (v-V)  p]   =   kV^T  + 


1-1$ 


(C35) 
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D.       SELF-CHECK    QUIZ    NO.     5 

1.  As    a    consequence    of    steady    flow 

a.  the    velocity    is    not    a    function    of   the    position    of 
the    fluid    particle. 

b.  density    is    not    a    function    of   the    velocity. 

c.  velocity    is    not    a    function    of    time. 

d.  density    is    constant. 

2.  What    term(s)    would    you   eliminate    from   Eq.    022    if  you 
are    told    that    the    fluid/flow   is 

a.Inviscid. 

b  .   Incompress  ib le . 

3.  If  for  a  fluid  element  it  is  found  that  there  exists 

a  linear  relation  between  the  rate  of  strain  and  stress 
the  fluid  is  said  to  be 

a  .  E  u  1  e  r  i  a  n  . 

b  .  Newton i  an . 

c.  Lagrangian. 

d.  Isotropic. 

4.  From  the  differential  form  of  the  Momentum  Equation  we 
derived  several  equations  associated  with  great 
scientists  and  mathematicians  of  the  past.   The  expres- 
sions known  as  Euler,  Bernoulli  and  Navi er-Stokes 
Equations  have  in  common  that  they  are: 

a.  Energy  Equations. 

b.  Continuity  Equations. 

c.  Moment  Equations. 

d.  Equations  of  Motion. 


auof^i 


q    'ap.LS    pueq-:;L|6La    aq:^   ^o    :tse[ 
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(q) 
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5.   Match  with  the  name  of  the  equations  (Column  1)  derived 
in  this  supplement  all  the  assumptions  or  constraints 
made  in  its  derivation  (Column  2). 


Col umn  1 

a.  Navier-Stokes  Equations 

b.  Euler  Equation 

c .  Bernoul 1 i  Equati  on 


Col umn    2 

a.  Inviscid    fluid. 

b.  One-dimensional     flow. 

c.  Isotropic    fluid. 

d.  Incompressible    flow. 

e.  Newtonian    fluid. 

f.  Steady    flow. 

g.  No    velocity    gradients 


^    'p    ' q    '6    uo    e      :o 

•6    JO    e       :  q 

•  9    ' p    *3       :  e 
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